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I. INTRODUCTION 


Production of vector boson pairs at the LHC is an important process for a variety of reasons. They 
include interesting checks of the Standard Model, searches for anomalous vector boson couplings 
and studies of a dominant background to Higgs boson production. Since the vector boson pair 
production cross-section is currently known to next-to-next-to-leading order (NNLO) in perturba¬ 
tive QCD [1, 2|, it is tempting to conclude that theoretical understanding of this process is already 
good enough and further improvements are unnecessary. Unfortunately, such a conclusion would 
be premature. 

Indeed, since production of vector boson pairs occurs predominantly through the annihilation of 
quark-antiquark pairs, the one-loop gluon fusion matrix element contributes to pp —» V\ V 2 cross- 
section at NNLO for the first time. The increase in the pp —>• V\ V 2 total cross-section caused by 
the gluon fusion contribution is not very large - it is about five percent in case of the W + W~ final 


state [2]. However, this gluon fusion contribution is somewhat larger than the theoretical uncertainty 
in the production cross section, as estimated in recent NNLO QCD computations jl, 2]. Therefore, 
for these uncertainty estimates to be valid, the current calculations of the gluon fusion contribution 
to pp —>• V\V 2 must be accurate to 30 — 50 percent. As we explain in the next paragraph, it is 
conceivable that this is not the case. In addition, the relevance of the gluon-fusion contribution 
strongly depends on the applied cuts and selection criteria. For example, typical cuts used by 
ATLAS and CMS in Higgs bosons searches increase the fraction of gg —» V 1 V 2 contribution to the 
background cross-section to 0(10%) [ 3 ]. With more aggressive cuts, this fraction can increase to 
an astounding 30%, as shown for instance in Ref. [4]. Good theoretical control of the gluon fusion 
process becomes very important in this case. 

Since the gluon fusion contribution may play an important role in four-lepton production at hadron 
colliders, it is important to realize that radiative corrections to gg —> V\ V 2 can naturally be large. 
Indeed, gg —>• V\ V 2 is a process where an initial state with large color charge annihilates into 
a colorless final state. NLO QCD radiative corrections computed for similar processes, such as 
99 H, gg ->• HH, gg -s> HZ and gg ->• 77 US turned out to be as large as 50 — 100 percent. 
Therefore, to fully control the gluon fusion contribution to vector boson pair production, NLO QCD 
corrections to gg —» V 1 V 2 need to be calculated. The first step in this direction is the calculation of 
the two-loop (next-to-leading order) scattering amplitude gg —>• V\ V 2 . Once the two-loop scattering 
amplitude becomes available, it will be possible to compute the gluon fusion contribution to vector 
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boson pair production cross-section through NLO in perturbative QCD by supplementing it with 
the real emission contribution gg —>• V\ V 2 + g. 


The goal of the present paper is to present a calculation of the gg —» V\ V 2 helicity amplitudes medi¬ 
ated by loops of massless quarks. We choose to work within the massless approximation because in 
this case all the relevant two-loop integrals are known {fj-12] 1 and because in many cases massless 
quarks provide significant fraction of the full amplitude. We note, however, that there are phe¬ 
nomenological situations where top quark loops are important contributors to gg —>• V\ V 2 fis| : our 
calculation is not applicable to those cases since for top quarks the massless approximation is obvi¬ 
ously invalid. We note in this regard that a proper treatment of massive quarks in loop-mediated 
processes is an important and difficult problem whose solution is not known at the moment. We 
therefore start with the massless quark contribution to the gluon fusion process gg —>• V 1 V 2 which, 
at the very least, will allow an informed estimate of the magnitude of NLO QCD corrections. 
Our result can then be supplemented with the contribution of the massive top quark loop once it 
becomes available. First steps in that direction where recently reported in Refs. [16]. 


The remainder of the paper is organized as follows. In Section Q7I we discuss the set-up of the 
calculation and the paranretrization of the gg —>• V\ V 2 amplitude. In Section IIIII we describe 
numerical implementation and checks of the calculated amplitudes and present some numerical 
results. We conclude in Section IIVI Finally, we note that many aspects of the calculation that we 
report here are similar to what we described in Ref. [I?!; nevertheless we believe that there are 
essential differences in the calculation of qq —>• V 1 V 2 [l2, 17] and gg —>• V] V 2 amplitudes at two loops 
to warrant a discussion of the latter in a separate publication. 


II. THE SET UP OF THE COMPUTATION 


We consider the process g(pi)g(p 2 ) —>• (Y]*(p:i) —>• Kp^)Kp &)){^2 (pA) KP7)Kps))- We work in 
the approximation where quarks of the first two generations are massless and quarks of the third 
generation are neglected. The CKM matrix is taken to be an identity matrix. We write the matrix 
element as 


M(A 91 ,A 32 ,A 5 ,A 7 ) = i(^j S a ^V^C^C^{\ b )el{\ 7 ) C Vl v, A^ 31 ,p* 32 ;P 3 ,p 4 ), 


(1) 


1 For earlier results in the case of equal-mass vector bosons, see EQ. 
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where a 1.2 are the color indices of the incoming gluons, gw = e/ sin 9w is the SU( 2 ) weak coupling, 
T>i = 1 /{pi — my. + imyTyJ is the V^-boson propagator, A sliS2 , A 5 , A 7 are helicities of the incoming 
gluons and outgoing leptons, Cfy C^y are helicity-dependent couplings of vector bosons to leptons, 
and £3 4 are matrix elements for leptonic decays of 1 \ and V 2 that we will specify shortly. The factor 
Cy 1 Vi involves sums of couplings of virtual fermions to gauge bosons. 


We now give these couplings explicitly. The couplings of electroweak vector bosons to external 
leptons read 


C^ R = -V2Q lS md w , c\f = 1 (t5±A t ), Cf w+ = Cf w - = S XL . ( 2 ) 

’ v 2 cos thy 

Here, V e = —1/2 + 2sin 2 dw, V u = 1/2, A e = —1/2, A u = 1/2 and Qi is the lepton electric charge 
in units of the positron charge. The couplings to virtual quarks, where two massless generations 
are included, are 


^ _ 20 sin 

°77 — g > 

2 sin 6w 


C'zz = 


{v* + vi + AI + AI) 


cos 


(3) 


Cz'y = — 


(VuQu + VdQd.) j C w +W - 


cos 0w 

where Q U) d are the electric charges of up and down quarks in units of the positron charge and 
V u = 1/2 - 4/3 sin 2 6 W ,A U = 1/2, V d = -1/2 + 2/3 sin 2 0 W . ; A d = -1/2. 


A word of caution is required here. As can be seen from Eq. (j3]) , the amplitude for the gluon 
fusion into a pair of vector bosons is proportional to sums of squares of vector and axial-vector 
couplings of these vector bosons to internal fermions. This means that diagrams with two vector 
currents and diagrams with two axial currents give identical contributions to the amplitude and 
that diagrams that involve one vector current and one axial current do not contribute at all. These 
two features have different origins. The first one is the result of the massless approximation and 
an ensuing relations between matrix elements of vector and axial currents. The second feature is 
a direct consequence of C'-parity conservation fl3 [ . that enforces the cancellation of these axial- 
vector terms when contributions of all diagrams are summed up. Therefore, the knowledge of the 
amplitude for vector-like couplings of gauge bosons to quarks is sufficient to reconstruct the gluon 
fusion amplitude for electroweak gauge bosons whose couplings to fermions are linear combinations 
of vector and axial-vector couplings. 


The scattering amplitude in Eq. CD stripped of all non-essential couplings reads 


^4 


A 

c 4 


— A c v 

— Sl/iuaf 3€i j A 1 e 2,A 2 e 3 e 4 > 


(4) 
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where 61,2 are polarization vectors for gluons with definite helicities. As we just explained, we 
need to compute the amplitude A for vector-like interactions of massive vector bosons. We need to 
decompose the amplitude A into invariant form factors that are independent of polarization vectors. 
To do this, we choose physical polarizations for gluons and vector bosons, 

ei • Pi ,2 = 0, e 2 ■ Pi ,2 = 0, e 3 -p 3 = 0, e 4 • p 4 = 0. (5) 

It follows from Eq.(|5]) that we use pi^ as a reference vector for polarization vectors £ 2 , 1 , respectively. 

To reduce the number of independent scalar products, we employ the Sudakov decomposition of 
the vector boson momenta, p 3 , 4 = a 3 , 4 pi + /3 3 , 4 p2 =tp_L- The coefficients 013,4 and /3 3 , 4 are expressed 
in terms of Mandelstam invariants 

mq — 74 „ ?? 7 q — t ml — t „ 777.4 — u 

a 3 = —2 -, (3 3 = —2 -, a 4 = —-, j3 4 = — - 

s s s s 


Since e 3 ■ p 3 = 0 and e 4 • p 4 = 0, we can choose £ 3,4 • p\ and £ 3,4 • P 2 as independent scalar products; 
£3 4 • p± are then given by linear combinations of £3 4 ■ p\, 2 - Given these constraints, the scattering 
amplitude is represented by the following expression 


A = T, (ei • e 2 ) (£3 • £4) + T 2 (£1 • £3) (e 2 • £4) + T 3 (£1 • £4) (f 2 • £3) 

+ (ei • £2) (T 4 (pi • £3) (pi ■ e 4 ) + T 5 (pi ■ £3) (p 2 • £4) + T 6 (p 2 • £3) (Pi ■ £4) + T 7 (p 2 ■ £3) (P2 ■ e 4 )) 

+ (ei • £3) (P± ■ £2) (Ts(pi ■ £4) + T g (p 2 ■ £4)) + (£1 • £4) (p± ■ £2) (Tw(pi ■ £3) + Tn (p2 ■ £3)) 

(7) 

+ (£2 • £3) (p± ■ £l) (TuiPl ■ £4) + T13 (p 2 ■ £4)) + (£2 • £4) (p± ■ ei) (T14 (pi • £3) + T15 (p 2 ■ £3)) 

+ (ei • P±) (£2 • P±) (Tnipi ■ £3) (pi ■ £4) + TisOi • £ 3 ) (p 2 ■ £4) + Ti 9 (p 2 • £3) (pi ■ £4) 

+T 20 (P2 • £3) (P2 ■ £4)) + (£3 • £4) (p± ■ ei) (p± ■ £2) Ti 6 . 

The form-factors T 1..20 are functions of Mandelstam kinematic variables s, t, u and invariant masses 
of the two vector bosons 7773 , 4 . 


The problem with the calculation of the scattering amplitude “as it is” is that carrying around 
polarization vectors for external gluons and vector bosons is extremely expensive for computations 
that employ integration-by-parts identities jl9, 2 f]], the main vehicle for multi-loop calculations. 
For this reason, we need a procedure that allows us to compute amplitudes but where no vectors, 
except for momenta of external particles, appear in the calculation. A suitable method employs 
projection operators. The idea is to compute the amplitude A in Eq. ([7]) with polarization vectors 
£14 substituted by linear combinations of external momenta pi, 2 , 3 , 4 - This procedure allows us to 
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compute the “projected” amplitude A in terms of form factors T 120 and also using integration-by- 
parts identities since polarization vectors completely disappear from the calculation. By choosing 
a sufficient number of independent “projection operators”, we produce a system of equations that 
we can solve for the form factors Tj 20 - 

However, there is a subtlety. The expression for the amplitude in Eq. ([T]) was written under the 
assumption that polarization vectors satisfy the transversality conditions Eq. (j5|) ; those conditions 
will be definitely violated if polarization vectors are replaced by linear combinations of external 
momenta. To get around this problem, we write the amplitude as 

A = A^r^ 1 V™ 1 V« ai vl Pl ei/u £ 2 ,, e 3oi e 4/3l , ( 8 ) 


where 


^nftu _ , P1P2+P1P2 , P 3 P 3 

7 12 — y ' 1 3 — y 1 

Pl * P 2 P 3 


spli'V _ M*' I P 4 P 4 

2 j 1 4 — y ' 2 ■ 

P 4 


(9) 


We then note that in Eq. ((HJ) we can use any vector to calculate the amplitude since projection 
operators Pi 2 , P 3 ,4 automatically ensure the transversality constraints 2 ■ E\ = 0 , 2 ■ 62 = 0 , 

P 3 • e 3 = 0 and ^ 4-64 = 0. We can also replace e\£ 2 e 3 e 4 with an arbitrary rank-four tensor since it 
will also be projected on the appropriate transverse space. 


We now list all projection operators that we use in the computation of gg —>• V 1 V 2 amplitude. To 
this end, we define the amplitude contracted with the projection operators 


Qn lA‘2A t 3M4 


12 7 12 7 3 7 4 


( 10 ) 


We also define a tensor that is a projector on the vector space that is orthogonal to the collision 
plane. It reads 


fM — XMPlP2P± nil 

V V U 12PlP2P± ’ V J 

where 5Zlplp $'2 = SvlfftXS pTpTpiPippP^^PPpa aIld 

«SSS 1 = det| s ^< I ‘;; 4 '‘ ) 1 |. ( 12 ) 

It is clear that contraction of tensor t fw with any linear combination of pi,P 2 and p± vanishes 
thanks to the antisymmetry of the determinant in Eq. (11211 . 


We define twenty projections of the amplitude A on linear combinations of Ti .,,20 by making different 
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choices of the “polarization vectors”. They are 


Gi = 0^^g piP2 g P3IM , 

q 2 = O^^g^g^, 

Gs = O^^g^g^, 

G 4 = pl 4 s~ 2 O p±p±pipi , 

G 5 =pi*s- 2 0 P ^ P1P2 , 

g 6 = p - 4 s - 2 o p ^ P2Pi , 

g 7 = p i 4 s~ 2 o p±p±P2P2 , 

G 8 = 4 pl 6 s- 2 O p ^ P3P % 3IM , 

G 9 = Apl e s- e O p ^ p H pilX2 , 

Gio = ^s- 3 O p ^ p H P2P3 , 

G u = 4 pI 6 s- 3 O p ^ p H P2P3 , 

G 1 2 = 8pl A s- 3 O p ^ p n flulli , 

Gy, = Apfs^O^^H^ 

G u = 8pi 4 s- 3 o^ p ^ p n PUP3 , 

G 15 = Apl e s- 3 O p ^ p H PliP3 , 

G 16 = 8pl 4 s- 3 O p ^ p H P2 , P4 , 

G17 = 4 p J 6 s~ 3 O p± P2P± jM4 , 

G 18 = Apl 6 s- e O pip2PiP2 t PllX2 , 

Gig = Apl 6 s- 6 0^ P2P2P H^ 2 , 

G2 0 = Apl 6 s- 6 O piP2P2P2 t PlP2 . 


( 13 ) 


In these equations, we used a simplified notation for the contraction of the tensor O with a vector a, 
_ Qui-a.-nn gj nce q 1 20 only depend on scalar products of external momenta and on 
scalar products of external momenta and the loop momenta, we can express G1..20 through known 
master integrals B, 10| by applying integration-by-parts identities Q 20 ]. At the same time, the 
form factors Tj .20 can he written as linear combinations of G1..20 in a straightforward way. 

To determine physical amplitude for the process g{p\ )g{p'i) —>• (V*(p 8 ) —> l{pb)KP&))(V*( pa) —> 
l(p7)l(p&)) we use spinor-helicity notations. Specihcally, for the incoming gluons we choose 

[27^1) M (2 7 ^1] p [1 7 ^2) „ <17^2] 

e 2 L ~ 


- _ 
e iL ~ 


\/ 2 [ 21 ]’ 


V / 2(21) ’ 


a/2 [12] ’ 


\/ 2 ( 12 )' 


( 14 ) 


Since a complex conjugation of the helicity amplitudes reverses all helicities, we only need to 
consider two, rather than four, cases of equal and unequal helicities. We choose L1L2 and L1R2 
as polarization states for gluons g, and g2, respectively. We leave the polarization vectors of the 
massive vector bosons unspecified at this point. To proceed further, it is convenient to write tensor 
products of gluon polarization vectors as follows ( see e.g. Ref. 4 j) 

( 12 ) 

ypi,nP2,w T pi,vP2,n - y/j,yp 1 ■ P2 -r ^ pupiP2 ), 

( 15 ) 




[12]s 
<lp±2] 


(Pl,pP2,v T Pl,vP2,n 9/i,vPl ' P2 T ^^lupipi) 
2 


[lp±2)p 2 ± .s 2 V 4 


P±.,/jJ>±vS 


+ sie piP2P±p p± l/ tpip2p±p. e pip2p±i' + (p v) 


The transverse momentum is introduced just before Eq. Q. We use Eq. ( 1151 ) to express the 
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amplitude A in Eq. fl7|) through nine independent Lorentz structures 


A\~l\2 N\i\2 


Fi lX2 (Pi ■ £4) (pi ■ e 3 ) + F 2 AlA2 (pi • e 4 ) (P2 ■ e 3 ) + E 3 lA2 (p 1 • e 3 ) (p 2 • e 4 ) 

+ -F]^ lA2 (p 2 • £ 4 ) ( P 2 • e 3 ) + -F]^ lA2 e 4 • e 3 + iifu/apPiP^Pl^ (^q iX2 Pi ' e 3 + F^p 2 ■ e 3 ) (16) 
+ i^ual3PlP2Pl4 ( F S lX2 Pl • £ 4 + ^9 : 1X2 P2 ’ € 4 ) 


In Eq. (ll6j) N\ 1 \ 2 are the normalization factors for left-left and left-right polarization cases 


Mll = 


( 12 ) 


Nlr = 


< 1 P± 2 ] 


[ 12 ]s ’ [lpx 2 )p^_s 2 ’ 


(17) 


and A2 9 are helicity-dependent form factors that are functions of the Mandelstam variables and 
the invariant masses of vector bosons. 

To account for transitions of vector bosons to final state leptons, their polarization vectors are 
replaced by matrix elements of vector and axial-vector currents. We therefore choose 


&, = <6i-m <&=<6i-m «jz, = wm =<siv*m. (is) 


We note that, although we need helicity amplitudes for all possible helicity combinations of leptons, 
it is sufficient to compute just one of them since other helicity amplitudes can be obtained using 
simple replacement rules. The amplitudes for left-handed polarization of both electroweak vector 
bosons read 

AlIl =AF AlA2 {(E AlA2 (15)[61] +E AiA2 (25)[62])(17)[81] 

+ (f 3 AiA2 (15)[61] +F AiA2 (25)[62])(27)[82] +2F AiA2 (57)[86] 

: (19) 

+ lA2 < 15 >[ 61 ] +E AiA2 (25)[62]) ((12)(78) [81] [82] + (17)(27) [21] [87]) 

- i(F AlA2 (17)[81] +F AiA2 (27)[82]) ((12)(56)[61][62] + (15)(25)[21][65]) }. 

Amplitudes for right-handed polarizations of the vector boson with momentum p% (p 4 ) are obtained 
from the above ones upon the replacement 5 -H- 6 (7 <->-8). Finally, all remaining helicity amplitudes 
can be obtain by replacing all angle brackets in spinor products with square brackets in Eg. (1191) 
and vice versa 

■A^rIr 2 = -4 3 l4L K*j) ^ K?]] ( 20 ) 


The F form factors that enter the amplitudes are expressed through either T or G form factors. 
This can be done in a straightforward way using Eqs. (171131151161) . Examples of corresponding 
relations are given in the Appendix. 
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Figure 1: Representative two-loop diagrams that describe production of vector boson pairs in gluon fusion. 


III. CALCULATION OF THE AMPLITUDE 


21] and process 


We apply the set-up described in the previous Section to the calculation of gluon-fusion amplitude. 
There are 93 non-vanishing two-loop diagrams that contribute to the gg —>• VV amplitude; some 
examples are shown in Fig. [T| We generate the relevant diagrams using QGRAF 
them with Maple and Form [22j. We compute the contribution of every diagram to the G and 
eventually F form factors. At this point, the result is expressed in terms of two-loop tensor integrals. 
These integrals can be classified in terms of six different topologies, three of which are planar 


and three are non-planar 

[ 


FIRE 


flO]. 


The tensor integrals are expressed through the master integrals 
computed in Refs. [9], |l£j], using integration-by-parts technology {0 . 2D], We employ the program 


23-f25| to achieve this. Combining contributions of different diagrams, we obtain the results 


for the eighteen form factors (nine for LL gluon helicity configuration and nine for LR gluon 
helicity configuration) that are required to describe all helicity amplitudes for gg —>• V) V 2 process. 
We note that, compared to the calculation of qq —>• V 1 V 2 amplitude, the case of gg —» U 1 V 2 requires 
more complicated reduction since tensor integrals of a higher rank appear. Nevertheless, FIRE can 
successfully deal with this challenge. 

As we already mentioned, the helicity amplitudes are expressed in terms of master integrals com¬ 
puted in Refs. 0, 10]. The analytic expressions for these master integrals involve various functions, 
including logarithms, polylogarithms of multiple ranks as well as generalized Goncharov polyloga- 


26] in the computer algebra 


rithms. To compute the latter, we use their numerical implementation 

program GiNaC ( 27 )]. We note that GiNaC can be called from both Mathematica and Fortran provid¬ 
ing multiple options for the numerical evaluation of the amplitude. 


The gg —» V\ V 2 amplitude appears for the first time at one loop; for this reason this amplitude 
is ultraviolet and infra-red finite. The two-loop gg —>• V\V 2 amplitude contains at most 0( 1/e 2 ) 
singularities, where e = (4 — d)/2 is the parameter of dimensional regularization. The divergences 
of the two-loop gg —>• Vi V 2 amplitude can be predicted in terms of the one-loop amplitude using 
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results of Ref. 


28] . The relation between one- and two-loop amplitudes becomes very simple if 


expressed through bare, rather than renormalized QCD coupling. It reads 


A 2 = -^e ei *A 1 + O(e 0 ), (21) 

where Ca = 3 is the QCD color factor and Ai t 2 are defined through the following expression for 
the amplitude 


A = y [s~ e Ai + a 0 s~ 2e A 2 + 0(a 2 s )] . (22) 

In Eq. (122 p we use ao = ai°*T(l + e)(47r) £ /(27r), where is the bare QCD coupling constant. 

A connection between divergences of the two-loop amplitude A 2 and the one-loop amplitude A\ 
given by Eq. (1211) is important for checks of the correctness of the calculation since the computation 
of A 2 proceeds without separation into divergent and convergent parts until the very end. 


We are now in position to present some numerical results for the gg -A V\ V 2 amplitude. To this 
end, we choose kinematics of an irreducible background to Higgs boson production and take the 
center-of-mass energy y/s to be the mass of the Higgs boson y/s = mu = 125 GeV. The invariant 
mass of the vector boson V\ is set to p 2 = with my/ = 80.419 GeV. The invariant mass of 
the second vector boson V 2 is set to 25 GeV. We take the vector boson scattering angle in the 
center-of-mass collision frame to be 7 r /3 radians. We also take decay angles of the lepton Z 5 in the 
rest frame of the boson V\ to be 65 = it/A and 925 = 7t/2 and decay angles of the lepton Z 7 in the 
rest frame of the boson V 2 to be #7 = 7t/6 and <^7 = n. The four-momenta of initial and final state 
particles are given by 


pi = (62.5,0,0,62.5), 
p 5 = (48.2561024468725, 
p 6 = (37.6127597971275, 
p-j = (19.5655688780000, 
p s = (19.5655688780000, 


p 2 = (62.5,0,0,-62.5), 

13.8697156788798, -28.4324101181205, 36.4400941989053), 
12.2010429705974, 28.4324101181205, -21.3881346746519), 
-19.2853793247386, 0, 3.29933778517879), 
-6.78537932473856, 0, -18.3512973094322). 


(23) 


Our results for leading and next-to-leading order helicity amplitudes at the kinematic point Eq. 
are shown in Table [0 The divergent terms of the next-to-leading order amplitude are compared 
to predictions based on Eq. (1211) and perfect agreement is found. We also compared the leading 


program MCFM 


order helicity amplitude with the results of previous computations [29 32], as implemented in the 


33], and found agreement. Finally, we note that we compared the numerical results 
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Helicity 

A\(e = 0) 

A^/A 1 (e = 0), 1/e 2 

A™/A 1 (e = 0), 1/e 

M( 2 )/Mi(e = 0), e° 

LLLL 

-5169.9932 + i 10017.414 

-3.0 

—9.45694415 — i 16.4895884 

42.4852911 — i 65.01495 

LRLL 

—6427.41534— i 2610.6160 

-3.0 

— 14.9422361 — i 9.2198662 

-6.774561932-f 71.66763 


Table I: Leading and next-to-leading order helicity amplitudes. Momenta of external particles are given in 
the main text of the paper. 



Figure 2: The ratio of finite parts of two- and one-loop helicity amplitudes R 2 = \A 2 (LL) + 

A 2 (LR)\/(\Ai(LL)\ + \Ai(LR)\) as a function of the vector-boson scattering angle. 

for helicity amplitudes reported in Table [J with the results of the independent calculation fl and 
found complete agreement. 


Finally, it is interesting to explore the numerical stability of gg —>• V 1 V 2 amplitudes that we com¬ 
puted in this paper. The numerical stability of such amplitudes is_ known to be a potentially 
sensitive issue as earlier one-loop studies showed, see e.g. Ref. 


32j 


35j. To study numerical sta¬ 


bility, we consider the same kinematic point as described above but we treat the vector boson 
scattering angle as a free parameter. We then compare the results of the double-precision imple¬ 
mentation of the form factors in a Fortran program with, effectively, arbitrary-precision calculation 
in Mathematica. We find that helicity amplitudes computed in these two different ways agree well 
for values of the scattering angle as small (large) as 6 = 1 ( 9 = 179) degrees. For such angles, the 
transverse momentum of a vector boson is just 0.5 GeV. We therefore conclude that our numerical 
implementation of helicity amplitudes is sufficiently stable to allow their use in realistic numerical 
calculations. As a further illustration of these numerical results, in Fig. [5] we show absolute values 
of helicity amplitudes as a function of the scattering angle. 
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IV. CONCLUSIONS 


In this paper, we computed the helicity amplitudes for the production of electroweak gauge bosons 
in gluon fusion gg —>• Vj Vj. The electroweak gauge bosons are allowed to have different masses and 
be off-shell; their decays to fermion pairs are taken into account explicitly. The helicity amplitudes 
for gg —>• V 1 V 2 are described by nine helicity-dependent form factors. We construct projection oper¬ 
ators to compute those form factors from Feynman diagrams, using integration-by-parts identities 
and the master integrals calculated by us previously. Analytic results for helicity amplitudes are 
implemented in a Fortran code that is available from the authors upon request. 


The results for the scattering amplitudes gg —>• W + W~ and gg —>• ZZ obtained in this paper open 
up an opportunity to compute the NLO QCD corrections to the production of a pair of electroweak 
vector boson in gluon fusion. Such calculations are interesting both at low and high center-of-mass 
collision energies. In the former case, they allow for a more accurate estimate of the irreducible 
background to Higgs production and the possible signal/background interference effects. In the 
latter case, precise predictions for gg —>• ZZ are important for improving prospects of constraining 
total Higgs boson decay width following Refs. 


studies 


36-38], as well as for more generic Higgs off-shell 


39]. The corresponding real emission corrections gg —>• ZZg have to be calcu 


case. Such calculations are clearly possible with the existing one-loop technology 
be done in the near future. 


ated in this 


40-43] and will 
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Appendix A: The relations between various form factors 


For the two possible helicity combinations, the invariant form factors, expressed through T form- 
factors, read 


f l l 


Fk L 


F 3 ll 


f[ l 


f 5 ll 


f 7 ll 


r LL 


f lr 


F LR 


[{ml - u ) (T12 + T 8 ) - (ml - t ) ( T 10 + T u ) - s(p 2 ± Ti 7 + 2 T 4 )] , 

[(ml — t) (T12 + Tg) — (ml — t ) (Tu + T15) — p 2 ± sTig — 2 .sTg + 2T2 + 2T3] , 


- [(ml - u ) (Ti 3 + Tg) - (m| - u) (T l0 + T u ) - p 2 ± sT 18 - 2 sT 5 + 2 T 2 + 2 T 3 ] 


2 [( m 3 — 0 (^13 + Tg) — (ml — it) (Tn + T 15 ) — . S (piT 20 + 2 T 7 )] , 

- ^ + 2Ti + T 2 + T 3 ) , F 6 ll = ( m 3 ~ ^ 2 (Ts ~ r2) - Tig + T u , 

i p^s 

(ml - t ) (T 3 - T 2 ) ll ( m l - *) ( T 3 - T 2 ) 

-till 1 15, -C 8 — --r 112 t g, 


p 2 ± S 


p±s 


(ml - u) (T 3 - T 2 ) 




+ T13 - T 9 ; 


^ [(771,4 “ t) (p]_s(T w + T i4 ) - 2 (ml - u ) (T 2 + T 3 )) + pl_s (p!_sTi 7 - (m§ - u) (T i2 + T 8 ))] , 
7 j ['(777,4 - 0 (p±s(Tu + T15) - 2 (777,3 “ 0 ( t 2 + T 3 )) - ((m3 - t) (T i2 + Tg) 


-p^.sTig - 2(T 2 + T 3 ))] , 

F 3 R = ^P±S (- (ml - u ) (T13 + Tg) + (mf - u) (T w + T u ) + p 2 ± sTi 8 ) 

+ (T 2 + T 3 ) (pl_s - (m§ - u) (ml - u)) , 

if* = ^(K - «) (pis(Tn + T 15 ) - 2 (m| - t) (T 2 + T 3 )) + pi.s(pisT 20 - (m| - t) (T 13 + Tg))), 
if* = ^pWiplTw -T 2 - T 3 ), if r = pls(T w + T U ) - (m 2 - u) (T 2 + T 3 ), 
if* = pls(T u + T 15 ) - (m| - t) (T 2 + T 3 ), if r = (ml - f) (T 2 + T 3 ) + pis(Ti 2 + T 8 ), 
if* = (ml - u) (T 2 + T 3 ) + p 2 ± s(T 13 + Tg). 


(Al) 


The square of the transverse momentum written through Mandelstam invariants reads p\ = — (tu — 
mlm 2 )/s. 

When the F form factors are expressed in terms of projections Gu,.. i2 o shown earlier, the results 
appear to be relatively simple. In particular, all spurious poles in d — 4 , present in the relations 
between T’s and G ’s cancel out in the relations between T’s and G' s. To give an example of these 
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relations, we show results for a few of the simplest form factors for the LL amplitude 

ll = -u){G l2 + <?ie) s{m 2 3 +pj)(Gi 3 + G 17 ) 

6 (d — 3)(m% — t) (d — 3) ( 771.3 — t) 

zriLL G 12 + Gi e + (t — m 2 )(Gi 3 + G\ 7 ) 

F 7 =- j— -, (A 2 ) 

ll _ ( m l - t)Gio - s(ml+p\)Gii Gi4 + {jn\ - t)Gi 5 
{d — 3)(m 2 — u) d — 3 

ll _ (tt - mpGu - Giq (zi - m\)G u + a(m| + p\)G^ 

9 d — 3 (d — 3)(m| — t) 
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